Low temperature electronic properties of Sr 2 Ru04 II: 
Superconductivity 
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The body centered tetragonal structure of S^RuCU gives 
rise to umklapp scattering enhanced inter-plane pair correla- 
tions in the d yz and d zx orbitals. Based on symmetry argu- 
ments, Hund's rule coupling, and a bosonized description of 
the in-plane electron correlations the superconducting order 
parameter is found to be a orbital-singlet spin-triplet with 
two spatial components. The spatial anisotropy is 7%. The 
different components of the order parameter give rise to two- 
dimensional gapless fluctuations. The phase transition is of 
third order. The temperature dependence of the pair den- 
sity, specific heat, NQR, Knight shift, and susceptibility are 
in agreement with experimental results. 



I. INTRODUCTION 

Sr2RuC>4 has quickly triggered a large interest because 
of its structural similarity to the cuprates and. the uncon- 
ventional superconductivity below T c *v 1 K.EJ The mate- 
rial is tetragonal at all temperatures a The three bands 
cutting the Fermi level with quasi two-dimensional Fermi 
surfacestrta can be mainbz.associated with the three ti g or- 
bitals of the Ru 4+ ionSjBp The transport is Fermi liquid 
like for T c < T < 30 Kl3 □ and strongly anisotropic along 
the c axis.EI The enhanced specific heat, magnetic sus- 
ceptibility and electronicLjmass indicate the presence of 
significant correlations .Ero£3 

The superconducting order parameter was proposed 
to jfL4¥£| p-wave symmetry analogous to superfluid 
3 He.E3t3 Tiiis-|Was supported by magnetic experimen- 
tal probes, but has not been proven unambiguously 
in spite of large effcpts-in Andreev spectroscc^ErlLS 
thermal conductivityj!3i23 and ac-susceptibilityc3 mea- 
surements. No indication for ferromagnetic, .correla- 
tions has been fopid.ETE3 The specific heat,c3~Ej ther- 
mal coiiductivityj^l and nuclear quadrupole resonance 
(NQR)EI are consistent with two-dimensional gapless fluc- 
tuations in the superconducting phase. For a more de- 
tailed overview see Refs. pp^-|34]. 

The present paper is part 11 of a series of three. 
In part I (Ref. |32| ) the quasi one-dimensionality of the 
kinetic energy of the d zx and d yz electrons has been 
used to derive a bosonized effective model. At inter- 
mediate coupling the interaction leads to a quasi one- 



dimensional model for the magnetic degrees of free- 
dom. The resulting spectrum of elementary excita- 
tions allows for the coherent description of the observed 
enhanced dynamical magnetic_susceptibilityc3 including 
the observed scale invariance.E£l Together with the two- 
dimensional correlations the normal state effective ma 
renormalization,EjO the specific heat coefficient ~t — 
and the static magnetic susceptibility enhancementsfilEr[L2 
can be described consistently. 

This work is devoted to the superconducting phase. 
The phase transition to the superconducting state in 
Si'2Ru04 is three dimensional and the specific heat has 
a shape that is similar to that of a, mean-field transition 
at the transition temperature.E3t!3 The low temperature 
specific-heat data are clearly inconsistent with the BCS 
theory.Effj This scenario suggests the combination of a 
mean-field induced three-dimensional transition (Sec. |J) 
with electronic correlations dominated by the quasi two- 
dimensional configuration (Sec. III). Similar approaches, 
have been successfully applied to coupled spin chains .Ell 
The model consistently can account for the experimental 
data concerning the temperature dependence of the pair 
density, specific heat, susceptibility and thermal conduc- 
tivity (Sec. |\]). A detailed comparison of the present 
approach with the p-wave picture is given (Sec. 

Part III (Ref. |d|) consistently explains the experimen- 
tally observed unconventional transitions under magnetic 
fields in terms of the model derived here. 



II. MEAN-FIELD APPROACH 

The generic model Hamiltonian is 

l,l' l,V,<7 

In this notation the electron creation and annihilation op- 
erators are c! and c, for orbital v with spin a on site 

I, ni V (J is the usual electronic density operator, t^, is 
the hopping matrix element. The intra-orbital Coulomb 
repulsion is larger than the inter-orbital repulsion, i.e., 

Ug£ =U >U 1 



U£££ and U > U 2 = Ktfi Thc 
Hamiltonian in Eq. (Q) is quasi two-dimensional.Ea 

The correlations in effective one-dimensional systems 
show power law behavior.o They are always more sin- 
gular than two-dimensional correlations which diverge at 



f 



most logarithmically.il Since the kinetic energy of the 
d zx and d yz electrons is quasi one-dimensieaal we expect 
their correlations to play a dominant role£3 

In order for the system to undergo a finite tempera- 
ture phase transition the coupling of the Ru02 planes 
is necessary.O The inter-plane d zv -d zv hopping has been 
estimated in Rcf. 32 to be t± » 0.02 eV and is an order 
of magnitude smaller than the in-plane hopping. The 
inter-plane hopping involving d xy orbitals is expected 
to be significantly smaller for geometry reasons.Q It fol- 
lows that the d zv -d zv subsystem should drive the tran- 
sition. Thjs_ conclusion is in contrast to most other 
approachesoO'C3c3 which assume the pair correlations 
of the d X y electrons to be dominant. 

The sensitivity of the superconducting instability on 
the dimensionality can be probed by experiments under 
pressure. Recent ultrasonic measurementsc£l suggest that 
the superconducting T c increases upon uniaxial pressure 
along the crystallographic c axis. This is reminiscent of 
an increase in the inter-plane coupling and is consistent 
with a pairing mechanism induced \m, inter-plane pair 
correlations as proposed by Baskaran£3 

On the other hand T c has been shown to decrease un- 
der hydrostatic pressure MP 7 ! This can be understood in 
terms of the decrease of the in-plane lattice parameter 
which in turn increases the in-plane hopping parameters 
which lowers the relative interaction strength and ren- 
ders the d zv -d zv subsystem more two-dimensional. It can 
thus be concluded that hydrostatic pressure lowers the 
in-plane correlations^ including those of Coopers pairs 
leading to a lower T c . Note that this assumption finds 
strong support in the decrease of the effective thermody- 
namic electronic masses underJxydrostatic pressure ob- 
served in dHvA measurements £3 



A. Normal phase instability 

In summary the discussion above suggests the presence 
of an inter-plane pair-correlation term between the pla- 
nar d zv -d zu subsystems driving the observed phase tran- 
sition. Such a term can be written in form of the lowest 
order coherent inter-plane pair-hopping term: 



E 



P 



(1,1') 



Z,/i.,£T Z',f,fT l,fJ, / ,<T r V 



(2) 



The orbital indices v, fi, v' , // S {x,y} are restricted by 
the Pauli principle for a' = a to p! ^ fi and v' ^ v. The 
second sum runs over nearest-neighbor sites where I and 
I are on neighboring planes. 

Note that such an inter-plane pair-correlation term 
may have different origins. Possible suggestions are 
Coulomb interactions screened on the scale of the inter- 
plane distance c/2 or electron-pihonon interaction where 
the phonon degrees of freedomEiJ have been integrated 
out. Since a discussion of the specific origin of the term 



would not be conclusive in the context of the present 
approach it is omitted. 

Instead, a phenomenological estimate for the pair- 
hopping amplitude can be obtained by requiring its value 
to be consistent with the single particle inter-plane hop- 
ping t±, i.e., V^' 1 * ~ t~\_/vp- The energy denominator is 
estimated by a typical in-plane energy scale. Here the ki- 
netic scale uf was chosen for the sake of definiteness. De- 
pending on the specific mechanism, an interaction energy 
scale U£'£, as the energy denominator is also conceivable. 
Since the system has been shown to be in the interme- 
diate couplingEl regime, i.e., vf ~ U£'£, , the numerical 
outcome is similar. Considering that t± w 0.02 eV0 
and vf = vf/cl r* 0.7 eVl23 one has t]_/vF « 6 K < «f 
suggesting that a mean-field decoupling of Eq. ^ is rea- 
sonable. 

Fourier transforming the Fermi operators in Eq. (|^) 
leads to the Hamiltonian 



Va 



/ j jy ^k,ii,o- C q — k,^' .a'^q— k' ,i>' ,a' C k' ,zact 
»<»)° k,k',q 



(3) 



with the effective Potentials 



t, a a c 

V q w — cos ~q x cos ~q y cos -q z 
vf z z z 
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— x, — y, — z 



(4) 



(5) 



the potential is extremal and attractive, i.e., V q ~ 
— t\/vF < 0. As a consequence of the body centered 
tetragonal lattice structure the are reciprocal lattice 
vectors. In other words, the body centered tetragonal 
lattice symmetry allows for an umklapp scattering driven 
pair instability. 

Within the configuration space of the d zx -d yz subsys- 
tem six possible pairing states can be found. There is 
an "equal-flavor" doublet which can be identified as a 
spin-singlet 



fb,±(<Z) = E ( C q-k,x,] C k,x,l ± C q-k,yA 1 



k, v ,l) ' ( 6 ) 



as well as a "mixed-flavor" doublet with a spin-singlet 
and a5 z = spin-triplet component 

Pl,± («) = E ( C 9-fc,*,T c k,y,i ± C q~k,x,l c fc, y , T ) ■ ( 7 ) 



The "mixed-flavor" doublet with spin-triplet components 
with magnetization quantum number \S Z \ = 1 is 



P2,±(q) — ^2 \ C q-k,x,] C fc,a,T ^ C q-k,x,l C k,y,\ 



(8) 



2 



The six components can be regroupedS as an orbital- 
triplet spin-singlet 



dependent weighing factor g\ of order one in the denom- 
inator of Eq. (|l2j) : 



P t }={p\ 



pt pt 
* o,-> r \.- 



and a orbital-singlet spin-triplet 

pst { pt pt pt 

1 t — I ■ r 2,+) ^2, — J 



(9) 



(10) 



(13) 



given here in vectorial representation.^ The pair op- 
erator components have a very similar structure as 
the p-wswe components initially discussed by Rice and 
Sigrist.Eil Notably, the P2,± components have the same 
spin-structure as the unitary state with E u .-symmetry 
proposed in the p-wave scenario for Sr 2 Ru0 4 Ej The de- 
grees of freedomJhxm a two-dimensional analog to those 
of the A phaseO'EJ of 3 He. The crucial difference is 
that all components P a ,± render even parity real-space 
wavefunctions because the antisymmetry of the total 
electronic wavefunction is established through either the 
spin-singlet or the orbital-singlet configuration. A closer 
account of the resulting symmetries is given separately 
in Ref. |§ 

Introducing the variables a E {0, 1, 2} and b E {+, — } 
and switching to time-dependent pair operators in the 
Heisenberg representation the bare pair susceptibilities 



The pair-hopping process in Eq. (g) with /i = //, i.e., 
a = in Eq. (|13|), requires the orbital /i to be unoc- 
cupied in the initial state which is very unlikely with 
rin ~ 1.3. Pair-hopping processes with /j, ^ // in Eq. 
(g), where a = 1 or a = 2 in Eq. (|l3|), are much less 
effected by the Pauli exclusion principle. Consequently 
9o ^ 9\ ~ 92 ~ 1 andr-the phase transition occurs in the 
mixed flavor channels E3 This conclusion is only altered if 

the in-plane equal-flavor-doublet correlations Xo b (Qi 1 ^) 
are much larger than those of the mixed-flavor sector for 
which there is no indication to be found. 

In the presence of Hund's rule coupling the spin-triplet 
energy expectation-values are always smaller than those 
of the spin-singlet £3 



»ti 



H IP: 



< \{Pl f \H\P 



(14) 



Consequently (P c 



= ( p c [-) = i p i,-) = °- Consistent 
with S = 1 moments on Ru 4+ impurities in Sr2lr04,t2l 
one finds mixed-orbital spin-triplet pairing or, equiva- 
lently, orbital-singlet spin-triplet pairing. The initial 
SU(2)<£>SU(2) symmetry of the two electron spins is bro- 
ken down to SO (3). The order parameter carries spin 
one aad explains the absence of a change-pf the Knight 
shifta and of the magnetic susceptibilitytj in the super- 
conducting phase. 

At this point it is important to notice that the es- 
timated value for Hund's rule couplingcj in Sr2Ru04 of 
Jh ~ 0.2 — 0.4jey is larger than the estimate for the spin- 
orbit coupling^] of A sa 0.1 eV. Consequently orbital- 
singlet pairing is possible even if the degeneracy of the 
The mean-field phase transition occurs at T — T c when d zx and d yz orbitals is lifted due to spin-orbit coupling be- 
cause the larger Hund's rule coupling over-compensates 
the effect. 



The random phase approximation (RPA) is equivalent 
in its static limit to the mean-field approximation. The 
RPA pair susceptibility is 



a , b 



the Stoner criterion x2(g if 0)| Tc = X^ fc J (0, 0)|r c 



1 + ^x2 



(12) 



-V~ 



V n is fulfilled. The estimated value of Vn ~ 6 



■q j - "0 10 1U11111CU - nicrauiiittKU »fl,nicu. ^ 

K is consistent with the low critical temperature of 
T c = 1.5 K. 

In order to make a quantitative prediction of T c the 
Xab (0' 0) have to be determined as a function of temper- 
ature within the effectively two-dimensional model based 
on Eq. ([!]) . This is not straightforward as becomes clear 
from the discussions in Ref. [32] and is left to future stud- 
ies. 

Mean-field approaches do not allow to control rigor- 
ously the locaLconstraint imposed by the Pauli excltu 
sion principled Since the filling factors of the orbitalscll 
are ~ 1.3 it is important in a mean-field approach 
to include effects of the Pauli exclusion principle at 
least in an approximate manner. Formally this can 
be achieved by introducing a phenomenological orbital- 



B. Comment on competing energies 

Since in the context of the superconducting cuprates 
it has become habitual to discuss the pair instability in 
terms of— the competition between kinetic and potential 
energiesEil it is useful here to clarify the situation in the 
present approach to Sr 2 Ru04. 

In the superconducting state of a BCS superconduc- 
tor the single particle excitations are gapped and hence 
the system looses kinetic energy in all spatial directions. 
The loss is over-compensated by the gain in condensation 
energy. The intuitive idea is that the system gains "po- 
tential energy" trough the attractive potential mediated 
by the phonons. 
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Here the situation is similar, the system looses kinetic 
energy both in-plane as well as out-of-plane. Iruthe light 
of the two-dimensionality of the band structured the loss 
of kinetic energy in the Ru02 planes is at least an or- 
der of magnitude larger then that along the c axis. The 
gain in pair hopping energy that one might anticipate 
through the inter-plane coupling—term Eq. (|^) is quite 
small — unlike the case discussedo in the superconduct- 
ing cuprates. 

Instead, the gain in condensation energy over- 
compensates the loss of kinetic energy as becomes ob- 
vious from the sine-Gordon model onto which part of the 
degrees of freedom in Sr2Ru04 are mapped in Sec. [II D: 



the ground state energy is lowered proportional to a given 
power of the gap even though the singleparticle excita- 
tions are gapped [c.f. Eq. and Ref. p3j. Or, viewed 
from the perspective of the collective condensation, the 
ground state energy is lowered because the single particle 
excitations are gapped. 

Moreover, the notion of potential and kinetic energy is 
obscured in the bosonized approach introduced in Ref. 

Notably in the one- 



32] and applied here in Sec. Ill 



dimensional limit with preserved SU(2) invariants the 
Bose fields and their dual fields are equivalentLX L3 and 
so are kinetic and potential energy. In this light the dis- 
crimination between the two appears obsolete and the 
energy gain is explained most naturally simply by the 
condensation of the electrons in the collective state be- 
low the gap. 



C. Superconducting phase 

The order parameter is a spin-one SO (3) triplet with 
components (Pf v ) ^ 0. Together with the global phase 
degree of freedom this implies that the to order parameter 
has SO(3)®SO(2) symmetry. 

These results are consistent with the implications from 
the bosonized microscopic model of the in-plane corre- 
lations discussed in Sec. III. There it will be shown 



that there are two degenerate saddle points with slightly 
different spatial anisotropies. This leads to a total 
SO(3)<£>SO(2)<£>SO(2) symmetry of the order parameter. 
The thermodynamic implications of the resulting gapless 
modes of the order parameter are discussed in detail in 
Sec. fill. 



The total Hamiltonian H+H m f usually permits a Landau 
expansi on of the free energy of the system. This is done 
in Sec. IIIF consistent with the bosonized microscopic 



model for the in-plane correlations to allow for a more 
quantitative analysis. 

An important critique of the orbital-singlet pairing 
mechanism is that in the absence of interaction electrons 
can only pair where the idealized, one-dimensional bands 
of the d zx and d yz orbitals cross in reciprocal space. The 
resulting pairing phase space is much too small to ac- 
count for the spec ific heat anomaly at the phase transi- 
tion [c.f. Sec. Jiy B |. At intermediate coupling, as relevant 
for Sr2Ru04.E^I interaction effects significantly enhance 
the pairing phase space. A close discussion is given in 
Ref. where it is shown that the interactions lead to a 



rather homogeneous gap function. See also Sec. IIIF 



D. 7 sheet: inter-band proximity effect 

The coupling of the different bands through the inter- 
band proximity effect and the resulting single transition 
temperature for all electronic-degrees of freedom has been 
discussed in earlier workJjOa In the present model the 
driving correlations are those of the d zx and d yz electrons 
as a result of the geometrically strongly enhanced iftter- 
plane coupling of with respect to the d xy electrons .ErE3 

To obtain a qualitative theoretical picture including 
the d xy orbitals one must consider the terms in the full 



Hamiltonian Eq. 



that arc left out of the d zr -d. 



Hamiltonian Eqs. (19) and (E3[) discussed in Sec. 



Ill 



H t=Yl l V C l 7 ,a C l',-y,a + U ^.7,1^,7,1- 

,cr I, a 

^ {Uj ni^^ni^^. + ni a , a ni^^) (17) 



1,1' a 



The hopping parameters where given in Ref. p! 



-0-39 eV, i™ 



t7',7 - = -0.27 eV, and t7'Z_- 

l 'tr 3 



The usual mean-field decoupling of the pair operators 



l,l+y ^ " ' """" "l.H-x+y 

0.11 eV. From the symmetry of the t^ g orbitalscl follows 
that in lowest order the interaction parameters E/q, fT 7 , 
and U2 are similar to those introduced in Ref. |3^ in the 
context of the d zx -d yz subsystem and are thus smaller 
but of the similar magnitude as vp . 

The inter-orbital terms ~ and ~ U] can formally 
be mean-field decoupled via 



pjtp- 



pjt ( pj 



pi {pi y-\(p s t )| 2 . (15) 



The mean-field decoupling Eq. (|15|) has to be applied 
to the inter-plane Hamiltonian Eq. (^|). The mean- field 
contribution H m { sa H± is 



2RePf(l)(Pt 



(16) 



+ ( C i,7,cr C J,iy,cr') A^,aA,v,a' 



( C l, 7 ,<7 C U,<T<)( C l,i,v C l,u,<r') ■ ( 18 ) 



The mean-field contributions (cj i7)CT Q )I/j0 -') couple the d xy 
electrons directly to the pair instability driven by the 
correlations in the d zx and d yz bands and induce pairs in 
the 7 band. There is only one phase transition. 
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Since, the inter-band proximity effect is usually 
stronger a single gap is assumed herein. It has been 
argued that Sr 2 Ru04 might represent a particular case 
where the inter-band pxoximity effect is suppressed^ 
leading to a double gapcll structure. No signature of a 
second gap has been observecLiii point contactEZl or ther- 
mal conductivity experimentstil though. 

In principle the inter-plane pair hopping term Eq. (^) 
must be extended to include the d xy or 7 orbitals, i.e., 



v, fx, v' , fi' S {x, y, 7} as well as the free energy Eq. (|58| 



Since from the orbital geometry g J 7 -C g 7 , 



9j,y < 



g x ^ y those contributions are small so they will only con- 
tribute small quantitative corrections. 



III. MICROSCOPIC IN-PLANE MODEL 



The microscopic model derived for the d zx -d yz sub- 
system in Ref. successfully has been applied to de- 
scribe the normal state properties in Sr2Ru04. Conse- 
quently it is reasonable to adapt the same approach for 
the in-plane correlations to describe the superconducting 
properties. To assure consistency in the notation of the 
present manuscript the elementary results are given in 
the following. Please refer to Ref. [32] for the details. 

As discussed in Sec. [n] the relevant correlations are 
those of the electrons in the d zx and d yz orbitals. The 
dominant hopping amplitudes are quasi one-dimensional 



and given by t — — if'f +y in Eq. (Q). The bands 

are linearized with Fermi velocity v F ss -\/3 to. To study 
the qualitative properties of the model with parameters 
relevant for Sr2Ru04 it proves useful to express the two 
spin and the two orbital degrees of freedom trough the 
charge (4> p ), spin (</> s ), flavor (4>i), and spin-flavor (c^ s f) 
Bose fields and their conjugate momenta (11^). The rep- 
resentation can be simplified by rotating the reference 
frame through x = (x + y) and y — -^(x — y) with 

r = (x,yY. The charge (magnetic) sector Hamiltonian 
including forward scattering terms is given by 



iJ c(m) = lim i 



' n f(sf) 



St) J 



+ K(m) [dx<t>p(s) + dy4>l( 
+ K(m) [<%^p(s) + dx(f>i( s f)] 2 } , (19) 



the spin (/i = s, v = y) and spin-flavor (/1 = sf , v = x) 
channels decouple. 



(22) 



It is then possible to approximate iJ m H s + H s f. The 
Luttinger liquid parameter K p and the velocity v p are 
effective parameters of the theory. In the case of SU(2) 
symmetry in the spin subspace the interaction term Eq. 
@ yields a rescaled -> K* = 1. 

It must be emphasized that the one-dimensional ide- 
alization of the magnetic subsystem albeit practical has 
to be used with caution. As discussed in detail in Ref. 
|32] in the Fermi liquid regime for T < 25 K ~ V m the 
magnetic correlations are two-dimensional but with cor- 
rections imposed by the closeness to one-dimensionality. 
Here the one-dimensional expressions will be applied to 
make use of results from the literature available for those 
systems but a discussion of the expected applicability to 
the realistic case including the two-dimensional coupling 



is given at the same time (Sees. HID and HI El). 

The back and umklapp scattering term in the 
bosonized Hamiltonian is 



Hint = 



Un 



d r cos v 47r0 s (r) cos v 47r</> s f (r) 



(27ra) 2 

+ - 1 .„ [d 2 r cos \y/i^<j) f (r) - 2\/2fc F yl 
(2iray J 

x ^ U\ cos VAnifis (r ) + [7 2 cos V&rrf'sf (^) ) 

+ J^ 7 cos [V4n<j> P (r) - 2V2k F x] 

x (ll 2 cos V47r0 s (r) + Ui cos VAn<j) s { (r)^ . 

(23) 

The limit a — > and L — > 00 is understood. Note that 
the symmetry of the superconducting saddle point (Sec. 
Ill B| ) deduced from the interaction term Eq. (|2^) is man- 
ifestly independent of corrections to the bosonized ap- 
proach since the energy scale is at least an order of mag- 
nitude larger than the correction terms, i.e., Uo.1,2 3> 

h,h,±,z > Vm M- Ref.|l. 

Equations ( |l9|) through ([23] ) are explicitly invariant un- 
der the symmetry operations of the tetragonal lattice. 



with 



Vc 



(m) 
(m) 



V F 
V F 



(-)Uo 
(-Wo 



[Ui + {-)U 2 ] 

Pi + {-)u 2 } 



(20) 
(21) 



e Hamiltonian of a crossed sliding Lut- 



Equation 
tingcr liquid. 

In the limit V m <C V m the magnetic correlations ex- 
hibit dominant onc-dimcnsionaLcontributions along the 
diagonals of the Ru0 2 planes.E3 On a mean-field level 



A. Bosonized pair operators 

At q — the pair operators defined in Eqs. (^|) - (||) 
are local, i.e., P a ,b{<l = 0) = Pa,b(l)- The correlations 
in the Po,± channel have been shown to be ne gligib le with 
respect to the mixed flavor channels in Sec. II A. In the 
continuum representation the pair operators can be given 
in terms of the Bose fields as introduced in Ref. B2L The 



abbreviations cos y / 7T(/) M (r) 



cos y/irQ^r) = cB^, sin ^/ttO ^(r) = s# M , as well as 
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4>i — \/2/7r /cf y and <j) p — 4> p — fcp £ are introduced 

for clarity. The pair operator components Eqs. (Q) and 
(H) then are 



P h+ {r) cx 



c9 s f c<t> { ccj) s + i s8 s { s(j)f s<f> a 



c9 s { c4> p c(j) si + i s9 s { s(j> p s<f> a t 



(24) 



Pi-(r) cx 



c9 s { s(j> { s(j> s + i s9 s [ c(j) { c<f) s 



+ c9 si s4> p s<j) s i + i s9 s f c4> p c(f> s{ (25) 



P 2 ,+ (r) cx 



no 



c9 s c(f) l ccf) s [ + i s9 s s0 f scj) s { 



+ c9 s c4> p c0 s + i s9 s scj) p sc, 



P 2 ,_(r) oc 



~2e~^e p 



TV (I 



c9 s s(f>f scj) s { + i s9 s ccjjf ccf> s { 



c9 s s4> p s(j) s + i s9 s c<j) p cq 



(26) 



(27) 



The limit a — > is understood. The Klein factors have 
not been plotted here for lucidity but are of crucial im- 
portance in the present model of multiple electron species 
when calculating correlation functions. They do not in- 
tervene with the symmetry arguments for the static ex- 
pectation values discussed below. 

A similar analysis of the spin density operators 



o.o 1 u=x,y 



(28) 



reveals that the out-of-plane component S z (l) ~ 

/(</> s ,(/) s f) depends on the Bose fields while the in-plane 

components S x (l) ~ S y (l) ~ f(9 s ,9 s {) depend on the 

dual fields. The a't , are the Pauli matrices. Con- 
ner 

sequently K* < K* f < 1 is consistent both with the 
observed magnetic easy plane configuration in Sr 2 RuOi 



(Sec. |VA] and Ref. |l|) and with (<!,")< Xy! (?> w ) ■ 



,( p ) 



X a b ' was defined in Eq. ([llj) . The bosonized model is con- 
sistent with the superconducting phase transition in the 
{P 2i+ ,P 2i _} sector. 

As discussed in Sec. II C of Ref. [32] the bosonized ap- 
proach neglects the exchange terms of Hund's rule cou- 
pling. This leads to an underestimation of the S z =JX 
spin-triplet component Pi +. The /iSR mcasurementsEj 
discussed in Sec. [VA and the upper critical fields (Ref. 
|l|) suggest that (P h+ ) « 0.8(P 2 ,+) 0.8(P 2 _). Since 

this is close to the isotropic limit Xi P +(Qi w ) ~ X^lily w ) 

will be assumed while x\ _ {<!■> w ) < Xi + (?> w ) an d conse- 
quently the spin-singlet component does not contribute, 
i.e., (Pi,-) = 0. The instability occurs in t he sp in-triplet 
channel consistent with the results in Sec. II A. 



B. Order parameter 

From Eqs. @, @ and @ follows that the finite 
value of the order parameter (P 2 .±), (Pi.+i-s^ implies 
(exp(— iy/7r9 p )) ^= 0. This is incompatible^ with finite 
expectation values of operators containing the dual field 
and thus terms ~ cos yfn <p p and ~ sin ypH 4> p are dis- 
carded. 

The interaction that led to the derivation of P m f in 
Eq. ( |l6| ) is attractive. A Tomonaga-Luttinger model 
with attractive interaction scales to strong coupling un- 
der a renormalizatiopLgxpup analysis and has a gapped 
excitation spectrum OO Its thermodynamic jpropjerties 
can be modeled via the sine-Gordon modelS Lj The 
free energy is lawjeced proportional to the square of the 
excitation gap.L3H The energy can be minimized by 
maximizing the gap. This suggests the minimization of 
the repulsive components of the interaction term H lnt 
through (cos y/^K(j) s ) = (cos */4m(f) s t ) — (cos -s/Itt 4>i ) = 
(cos v4tt 4> p ) = 0. It is then consistent to assume that 
these operators are irrelevant and consequently the con- 
tributions ~ cos^/nips, ~ cos v^^sf; ~ cos\/7r0f, and 
~ cos y/n^p m the operators P 2 ,± and Pi,+ can be ne- 
glected. [—. 

One then obtains a gap in all channels fi £ {c, s, sf}o 
with possible finite expectation values for (siny/w<p s f), 
(sin a/7t# s ), (cos tJtt9 s ), and (sin y/wfif) as will now be dis- 
cussed in detail. 

The system is invariant under a global phase shift 
in the charge sector. The fluctuations of the charge 
phase are Referred to as the Anderson-Goldstone collec- 
tive mode.E3 They can be parameterized by introducing 
the unit vector 



tin = 



-i sui^/tt9p 
cos ^/tt8 p 



(29) 



The Anderson-Goldstone collective mode results from a 
broken Gauge symmetry and in the presence of long 
range Coulomb interaction acquires a gap of the order 
of the plasma .frequency through the Andcrson-Higgs 
mechanism.OuJ The large gap allows for neglecting the 
charge phase fluctuations in the discussions of the low 
energy excitations that follow (Sec. HIE). 

Correspondingly, fluctuations of P 2 ,+ and P 2 ,_ can be 
parameterized as phase fluctuations in the spin channel 
via 9 S which takes the role of a azimuthal angle of the 
magnetic moment of the Cooper pairs. The third spin- 
triplet component Pi,+ can be included via a polar angle 



— i sin \Jtt9 s sin 9 Z \ 
cos y/n9 s sin 9 Z 
cos 9, 



(30) 



The SO (3) vector O s descri bes t he three spin-triplet com- 
ponents introduced in Sec. II A. 



G 



The order parameter expectation values 
(Pz,±), ( p i,+) ~ ( sm V^f) = (sin(V^>f - V2y)) break 
the invariance of the model under reflection y — ► —y . 
This in an artifact of the mean-field decoupling Eq. Jig ) 
which yields only one of two equivalent saddle points. In 
the absence of external fields, temperature gradients, or 
strain the symmetry must baxestored by including an- 
other two-component vectorErLJ 



n f = 



sin(77T0 f -y/2y)\ 
sin(y / 7r0f — \/2x) J 



(31) 



Equation (^) includes a p-wave component of the pair 
operator. This becomes apparent in the rewritten form 



cos V2y sin y/n^i 
cos \/2x sin y/itfa 



sin V2y cos y / 7r</>f 
sin y/2x cos y/n<f>f 



(32) 



Since the component f2f jJC is oscillatory along y the pair 
correlations are stronger along x and vice versa. The 
existence of such two order parameter components that 
couple differently to magnetic fields through a slight spa- 
tial anisotrqpy. has.been implied experimentally, via upper 
critical field^Oo and thermal c ondu ctivityEj measure- 
ments as discussed closer in Sec. IV E . The fluctuations 
of flf account for the observ ed str ain anomalies^ in the 
superconducting phase (Sec. HIE). 
Introducing the normalized average 



n^ = (n M )/|(n M )| 



(33) 



\flp\ = 1 is satisfied. Denoting the order parameter 
amplitude 

P=|(P}| = -|(sin^<M| |<n„>| I (a) I |(«f)| (34) 

7T 



the mean-field Hamiltonian Eq. (16) becomes in its con- 
tinuum representation 



H 



mi 



(35) 



— sin V^r^sf (si* p O p j (fl* s fl s ^ (fl* { flf 



Together with the in- plane contributions Eqs. (|19|) and 
( p2| ) the Hamiltonian in the superconducting phase is 
given by 



H sc — H c + H s + H s f + H, 



Imf ■ 



(36) 



C. Magnetic and temperature fields 



Magnetic fields need to be included in the microscopic 
model by introducing a vector potential in the momen- 
tum operator.L3 A detailed study within the model pro- 
posed here is rather involved and is left for future stud- 
ies. One can use the notion of spin-charge separation 



thought to discuss some qualitative implications of an 
applied magnetic field. The coupling of an applied mag- 
netic field h to the magnetic degrees of freedom of the 
order parameter then is simply given by 



H h = - 



d 2 r h fL 



(37) 



This term further reduces the symmetry in the spin chan- 
nel, e.g., fi s = (y/2 , \/2 , 0)t if the quantization axis 
is chosen along the applied field. The ground state order 
parameter then carries spin one along the direction of the 
magnetic field. 

The study of the coupling of an external magnetic field 
to the anisotropic components of the flavor vector flf is 
more involved. Related two-component superconductors, 
have been studied via a Landau Ginzburg analysiaZS'EZI 
albeit based on p-wave symmetries wi thout the notion of 
gapless fluctuations addressed in Sec. HIE . The general 
physical result though is that the magnetic field couples 
strongest to the order parameter component with dom- 
inant superconducting correlations perpendicular to the 
field. If the applied field has a component in the x-y 
plane, the coupling to the components of Of will be in- 
equivalent. A field along x will couple stronger to the 
order parameter component Slf tV since it has dominant 
pair correlations along y. 

The application of a temperature gradient in the x-y 
plane breaks the reflection symmetry of the lattice and 
yields inequivalent components of flf. Thepcder param- 
eter symmetry is reduced to a two-fold axisEil perpendic- 
ular to the x-y plane as opposed to a four-fold axis in 
the absence of the temperature gradient. Similar effects, 
are expected through strainEJ and at surfaces or grair£2l 
boundaries. 



D. Mapping onto sine-Gordon models 

The mean-field contribution Eq. ( |35| ) couples all the 
channels fi € {p, s,f, sf}. Since the lowest bound state 
and thus the gap energy is correctly reproduced by a per- 
turbativa-treatment even in one-dimensional interacting 
systemsEj one can attempt a mean-field like decoupling 



valid at least for small values of the fields 



To this end 



all terms ~ sin y/ir(f>p and ~ sin y/nOp are transformed to 
cos terms via the sliding transformations <f>p — > <j)p — ? 
and Op — ► Op — j. Expanding the cos terms to second 
order in the fields decouples H m f je ff in the different chan- 
nels [i. Rewriting the resulting contributions H m f jfi in 
terms of cos and reversing the sliding transformations 
Op — ► Op -\- ^ where necessary yields the desired result 
(up to a constant). 

In order to obtain the Lagrangian for the Hamiltoni- 
ans in Eq. (^2|) Euclidean time dependence of the fields 
is introduced, the Lagrange transformation Hp (r, r) = 
■^4>p(r,T) is performed, and the fields are rescaled as 
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The effective, one-dimensional ac- 



tion for the spin-flavor channel is obtained as 

L 

Srf =^rJ d yJ dT l#£ + (%^f ) 2 + FAf Bf cos B af ^ ] 



(38) 



in the absence of a gap. It is thus permissible to use 
the theoretical results given in the literature for the one- 
dimensional actions Eqs. (|3|) and ( |3^ ) and anticipating 
corrections from the two-dimensional character of Eq. 
( pfl| ) and the two-dimensional corrections^ to Eqs. (138 
and (H). 



with B s [ — \/i^K* f and inverse temperature (3 = T^ 1 . 
Equation ( p8[ ) is the quantum sine-Gordon action with 
the interaction parameter — v 8 ^° a 



_ 8Vo 
v^-ira 2 ' 

The transformation to the dual fields II,, 



-dxOft and 



rescaling as i9 M = (K*) 1 / 2 8^ yields the effective, one- 
dimensional spin-channel action 



v B L 





dx I dr 



i'i • i,M_: J • VMJIJK (39) 



-L 







for correspondingly adapted vectors Q s and Q' s . The 
action in Eq. ( |39| ) is invariant under a rotation of f2 s . 
Transversal fluctuations of £1^ with respect to Cl s give 
rise to Goldstone mode s whic h can be treated separately 
as will be done in Sec. HIE and Appendix A. To study 



thermodynamic properties as implied by the sine-Gordon 
action one can consequently use the simplifiecLfixpression 
with U' s n' s ->■ cosB s tf s with B s = ^/^/K*M 

For T — > and L — ► oo the sine- Gordon actions Eqs. 
(|38|) and (J3^) are equivalent to the charg£.jzjero-sector 



of the massive Thirring model if < 87r.EE For both 
channels to be in the gapped regime requires K* { < 8 and 

k;>1/ 

the conditions on K* are satisfied 
For T 



have a gap ofEfl 



Unless K* < 1 is very strongly renormalized 
re satisfied. 
00 the actions Eqs. (£j 



and 



A M ~ [P a 2 M M ] ' 



1VC47, 



(40) 



to the lowest bound state above the ground state. 

Similarly the action for the charge channel can be de- 
rived. 



E. Fluctuations 

The description of Goldstone modes in symmetry bro- 
ken phases via non-linear sigm a mod els is a standard 
problem in text book literature.E2ElllI3 The description 
of the charge phase fluctuations including the Anderson- 
Higgs mechanism leading to the charge-phase gap has 
been derived by Schulzc3 for the case of the attractive 
two-dimensional Hubbard model. Goldstone modes of 
the magnetically ordered phase of the repulsive two- 
dimensional Hubbard_model have been discussed by 
SchulzHH and WeneH et al. Most importantly, the 
present app roac h is a two-dimensional analog of super- 
fluid 3 He-AE3C3 The conclusive result is that for temper- 
atures sufficiently smaller than the electronic excitation 
gap T«A and at sufficiently small energies uj <§C A the 
gapless fluctuations of the internal degrees of freedom of 
the order parameter associated with continuously broken 
symmetries are well described by non-linear sigma mod- 
els (Appendix A). 

In the present approach the gapless Goldstone modes 
correspond to the broken local spin rotational invariance 
described by the SO (3) vector f2 s and the degenerate 
flavor saddle points described by the SO(2) vector fif 
with local constraints |J7 M (r,r)| 2 = 1. 

At sufficiently small temperatures and energy scales 

T, u> -c A the stiffness of the non-linear sigma 

model describing the angular fluctuations has been 

determinedCirE 2 ] to scale with the square of the order 

— 2 

parameter amplitude ~ P . The analogy to mag- 
netic excitations in the ordered Hubbard model suggest 



correctional, that depend on the specific correlations 
of the system. The resulting nonJinear sigma model in 
Euclidean space (Appendix A) isc3 



1 



S = \ J d z r ^ dr <j v F (<j> 2 p + 4>f) 



V C [Chc(j)p + dy(j)f} 2 + V C [8y(j) p + 8x(j){] 2 

. ^ 4V& 



(41) 



The action Eq. @ can be transformed to both dual 
fields P and 0t but a mixed representation in 9 P and fa 
as desirable is not possible. 

No detailed theoretical results are known for the two- 
dimensional two-component sine-Gordon model defined 
by Eq. (fid]). In the gapped phase relevant here the dif- 
ference between the two-dimensional correlations and the 



one-dimensional ones are less severe (Sec. Ill El) then 



Sn = P 2 



d 2 r 



-L 



drJ^A 

IS.fl 



W 2 

(av„) 2 



the stiffness parameters are 

d 2 



^=fl3(^xg (9,0) 



with the excitation velocities 
1 



q=0 



V 2 r) 2 

M 2 A v dcu 2X2 ' ±[ ' ' 



> (42) 



(43) 



(44) 



Here the magnetic channel is treated isotropically. The 
pair correlation functions depend on Qx,<ly ~ Qx ± %■ 
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Consequently A x — A y reflecting the symmetry of the 
lattice. The stiffness perpendicular to the xy plane is 
negligible since the out-of-plane correlations are two org 
ders of magnitude smaller as indicated by the resistivityQ 
thermal conductivity,! 2 ^ and critical magnetic fielcS mea- 
surements. Within the mean-field approach A z = 0. 
Note that Eqs. (|4|) through ( p| ) can also be obtained 

by replacing X 2 — > P in the amplitude fluctuation ex- 
pression Eq. in Appendix B. 

The crucial difference of the present model with re- 
spect to BCS theory is the presence of gapless modes 
in the ordered phase. Unlike in BCS the excitations 
of the internal degrees of freedom of the order param- 
eter, namely the SO (3) spin-triplet components and the 
SO (2) flavor components do not necessarily acquirp-a gap 
(Appendix .B2). The analysis of a related modeE3 and 
the analogycJtlj to 3 He-A suggests that the magnetic es 



for 



mBEl 



ss and have linear dispersions .Ell 
hold for the flavor degrees of 



citations remain gaj 
Similar aijgumentsE-! 
freedomES 

Note that even in the case of a large gap VIa > A 
of the magnetic fluctuations O s which could in princi- 
ple induced by the spin-flavor coupling terms in Eq. ( [23] ) 
the physical implications of the present approach remain 
unaltered because of the ungapped SO (2) flavor modeE3 
Of. The degeneracy of the components of Of are prereq- 
uisite for the presence of a gapless flavor mode. At this 
point it should be stressed that the tetragonal structure 
of Sr2Ru04 has been found to be very stable (c.f. Ref.g). 
Since Qf, x and Qf^y are related via a symmetry operation 
under which the tetragonal point group of the lattice is 
invariant it can be concluded that the degeneracy of the 
two components is also very stable, at least in the bulk 
material. For a discussion see Sec. @). 

Experimental evidence of the existence of the flavor 
mode stems from ultrasonic measurements in Ref. ^9] that 
probe the strain dependence of the superconductivity in 
Sr2Ru04. The in- plane anisotropy of the components 
Q{ ;X and Of (J/ leads to a coupling of the superconduct- 
ing order parameter to the strain and thus accounts for 
the observed effects. The coupling to out-of-plane shear 
strain components is much weaker. The gapless exci- 
tations yield a reduction of the elastic constants in the 
superconducting phase. The strain breaks the symme- 
try between the two components which accounts for the 
anomalies observed near the phase transition. Also, the 
observed anomalous absorption in cyclotron resonance 
experimental might be attributed to the coupling to the 
Of gapless mode. 



1. Estimates for the pair correlation function 

On the mean-field level one expects that the pair cor- 
relation functions that determine the the stiffness and 
excitation velocities through Eqs. (p3[) and (Eh have the 



XbcsW» w ) 



~A/(q) 



T? a 2 — A 2 ' 



(45) 



where A is the single particle energy gap, v c s is the ef- 
fective velocity of the collective excitations associated 
with the pair operators, and f(q) is a function that 
weakly depends on q. Note that in the static limit and 
with a pair correlation length £p related to the super- 
conducting gapEj as fp oc v e s/A the pair correlation 
function has Lorentzian form x^ p ' (q — > 0, ui — 0) ~ 
^p 1 (|aq| 2 + Cp 2 ) -1 ■ Equation fl45| ) is thus valid on a 
rather general basis. 

Since one expects that the quasi one-dimensional cor- 
relations exhibit the most singular behavior, and since 
the normal phase magnetic properties of the system are 
dominated by quasi one-dimensional correlations,!! 2 ] it is 
necessary to investigate the possible impact of the one- 
dimensional subsystem on the stiffness and excitation ve- 
locities through Eqs. ([43|) and Q). 

First of all it must be stated that thq-spale invariance of 
the magnetic correlations breaks dowr£3 in the Fermi liq- 
uid regime for T < 25 K indicating that the system turns 
two-dimensional at low temperatures S3 Furthermore, the 
pair correlation functions factorize into spin, spin-flavor, 
charge, and flavor contributions in the bosonized repre- 
sentation [c.f., Eqs. ( p4| ) through (|27|)]. In the super- 
conducting phase the correlation functions have to be 
evaluated with respect to the actions Eqs. (38), (|3S|), 
and (^l|) where all channels are gapped. Consequently in 
the superconducting phase the divergent correlations of 
the quasi one-dimensional magnetic subsystem are sup- 

pressed and \\ ±(9> C4 i s analytic in q and u>. The con- 
tributions from the charge and flavor channels are two 
dimensional anyway [Eq. (f4l|)]. 

To make these reflections more transparent consider 
that for T > T c the spin and spin-flavor channel cor- 
relation functions are approximately determined by the 
Hamiltonian Eq. (H|) . They have the standard functional 
dependence in conformal field theory.EilO Density and 
pair correlation functions are connected via substituting 
|2fcF — qjj\ — > qu with Fermi wavenumber fcci-and an ap- 
propriate adoption of the scaling dimension.o Dimerized 
spin chains are described by the action Eq. (38) witho 
B 2 . = 2ir and the density correlation functions have been 
given in Ref. [h| The_value of B 2 = 2tt is not a singular 
point in the analysisE^I at T = and similar behavior 
for other values of B 2 ^ is expected, especially at finite 
temperatures. This leads to the estimates 



A cos 2 (ag T7 /2) 
-v 2 cS sin 2 (aqtr) - A 2 



(46) 



where A is the energy gap to the lowest bound state. 
Comparing with Eq. (45) one consequently finds that 



in the gapped phase Xid(?k> w ) = Xbcs^"' w ) witn 







-A/(q-) = A cos 2 (ag F /2). Moreover, Xi D 0) is an- 
alytic in qjj = 0. This result is not surprising in the light 
of the relevant three-dimensional couplings that have to 
be present in order to drive apx quasi one-dimensional 
system into the ordered phase. c£l 

In Appendix Bl it is shown that the analytical prop- 

(p) 

erties of the pair correlation functions \\ ± are related 

( P) A 1 

to those of Xbcs [ c -f- E( l- Wl- Given the analyticity of 
the one-dimensional contribution Eq. ([l6|), the coupling 
termSr-that render the magnetic system two-dimensional 
belowE3 T < 25 K, and the two-dimensionality of the 
charge contribution as shown in Eq. (|4l| ) it can be con- 
cluded that the A v as defined in Eq. ( [43| ) are dominantly 
isotropic in the q x -q y plane for T <C v c s and q v <C a -1 
up to quadratic order in q v . 

From the geometry of the upper critical magnetic fields 
discussed in Ref. [IT] it can be concluded that there is an 
in-plane anisotropy in the order parameter of ~ 7% which 
should be reflected also in X2 ± ■ Most importantly, the 
non-analyticity along the diagonals of the basal planeEJ 

suggests the presence of a small term ~ Ife^yl in the 
(p) 

expansion of X2 ±- Considering the small energy scale 
set by the anisotropic component of the upper critical 
magnetic field H c , a (T -»• 0) ps 0.05 T - 0.04 K the effect 
should be negligible here for T > 0.1 K. The expected 
significant next higher order anisotropic correction terms 



are ~ q 2 x q 2 



2. Effective Saddle point 

The amplitude fluctuations of the order parameter be- 
come relevant on energies scales of the order of the ex- 
citation gap lo ~ A (or, equivalently, for wavevectors 
a\q u \ ~ A/tv) and qualitatively change the behavior of 
the gapless modes at these energy scales (Appendix B). 
The description of the latter via the non-linear sigma 
models in Eq. ([12]) is consequently restricted to small 
temperatures, small energies, and large wavelengths. Eil 

As the gap decreases while the system approaches the 
phase transition thermal fluctuations become sufficiently 
large to excite the amplitude mode which in turn cou- 
ples back to the angular Goldstone modes. Amplitude 
fluctuations at energy scales of the order of or larger 
than the gap can be treated outside the critical region 
by the standard approach^ where Gaussian fluctuations 
are integrated out. In Appendix B2 it is shown that the 
Gaussian amplitude fluctuations can be included in the 
description of the f2 M modes in Fourier space in form of 
a cutoff and a fluctuation dependent prefactor: 



rtamp 



(3N 
4^2 



(X 2 



0<w„<A 

£ 



A v d 2 [aq] 



[aqu] 



(47) 



The number of in-plane Ru ions is N = L 2 /a 2 and the 
prefactor (X 2 ) is the variance of the amplitude fluctua- 
tions and parameterizes the effect of the amplitude fluc- 
tuations on the action of the Goldstone modes. Follow- 
ing the analysis in Appendix B one can approximate 

(3(X 2 ) « Xbcs (0.0) = i [ E q- ©]■ Near the P hase 
transition, where A — > 0, the prefactor diverges jaihich 
is reminiscent of the diverging correlation lengtho For 
T —> T c one has S™ p > Sq and the impact of the fluc- 
tuations of the internal degrees of freedom of the order 
parameter is determined by Eq. (f|7j). 

A dimensional analysis suggests that the integrals in 

Eq. ([It]) are proportional to A 5 . Using A v w ^Mjt [Eq. 

(fjj)], and v v w v c g [Eq. (|S5l)1 the contribution to the ac- 
tion from the amplitude-Goldstone mode coupling term 
can thus be estimated to be 



rtamp 



A 3 + 0(A 5 ) 

47T 3 V 2 



(48) 



efi 



amp 



The approximation is justified when the stiffness p 

f3{X 2 )A u w T^f is large enough, i.e., in the limit A 0. 
In this approximation all channels fi yield the same con- 
tribution. The dimensionless prefactor Sq ~ 1 is phe- 
nomenological since the exact numerical prefactor is not 
known within the approximations made. 

Equation ( [48| ) is an effective saddle point contribution 
in which the fluctuations of the internal degrees of free- 
dom of the order parameter have effectively been inte- 
grated out. The approach is only applicable outside the 
critical region T < T c — 0.02 K( Appendix B3). 

The contribution from Eq. |48| is of third order in the 
gap A 3 near the phase transition. This an important 
consequence of the two dimensionality of the problem. In 
one or three dimensions the contribution is A 2 or A 4 , re- 
spectively, and consequently only renormalizes the usual 
parametersO of a Landau expansion. 



F. Free energy 

The general expression of the total free energy con- 
sists of the contributions from the spin, spin-flavor, and 
charge channels and the mean-field contribution, i.e., 
F = ^2^Fi^ + FmS, with 



= -Tin Jdip^e s " for ^e{sf,c}, 



(49) 



F^ = -Tin Id&n e~ s » for fi E {s}, (50) 



where du 



6f. The actions are given by Eqs. 



(||), ©, and @ with (P) = in the normal phase. 
The mean-field contribution is 



F mf = -V NP -Tin 



n 



q,n(2n) (47r) 



-Sn 



(51) 
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The second term accounts for the gapless fluctuations of 
the i ntern al degrees of freedom of the order parameter 
(Sec. HIE and Appendix). 

The order parameter is obtained by minimizing the 
free energy via the variation (6F)/(6P) = 0. At finite 
temperatures the analysis of the expectation values pf the 
sine- Gordon and related models is quite involved. L3 An 
estimate of the symmetry of the full Eliashberg equation 
is discussed phenomenologically in Ref. [54|. Based on the 
result that interaction effects can render the gap function 
rather -homogeneous in agreement with the experimental 
results^ (Sec. |lVE| ) one can derive a Landau theory for 
the free energy near the phase transition in the ordered 
phase. 

To this end we omit at the moment the contribu- 
tion from the fluctuation action Sq, assume an isotropic 
order parameter ,00 and consider the resulting self- 
consistency equation which reads 



(52) 



L 2 V 



To obtain an estimate for Eq. 
equivalence of the actions Eq 




we make use of the 
and (I39h to those 
The self- 



of spin chains with alternating interactions 
consistency Eq. (^2|) is thus equivalent to that studied 
in the context of the spin-Peierls transition in CuGe03 
studied in Ref. |9(| The results show that a one- 
dimensional system that undergoes a phase transition 
to long range order — induced by an effective three- 
dimensional coupling that can be treated mean-field 
like — can be described near the phase transition by the 
standard mean-field free energy functionaE3 



F„ + R 



mf 



T~T C 



z — ^ a. 



NV Q 



T 



NbP 



(53) 



which is the value for free n 
Eq. ([io| ) then yields a gap oil 



issive fermions.O Near T, 



A fa 8 V P . 



(57) 



As discussed in Sec. [HE 2 the the gapless fluctuation 
contribution at the phase transition can be described via 
the effective saddle point term Eq. ([f8|), i.e., Sn — ► S^ 119 
in Eq. (|5l]). Together with Eq. ( p33] ) the Landau expan- 
sion of the free energy Fa near the phase transition then 
is 



Fa 

N 



N 



V tP 



3„2 



6P 



. (58) 



J off 



where t = 1 — T/T c <C 1 is the reduced temperature. The 
approach is valid outside the critical region, i.e., for \t\ > 
0.01 as estimated in Appendix B3. Since the expansion 
is in P > and not in the complex order parameter the 
phase transition described by Eq. iEa) is of third order in 
the sense of Ehrenfest's definition.E2JGauge invariance is 
preserved since P = \e l< ^ G (P)\ for arbitrary gauge fields 

4>G- 



IV. APPLICATION TO Sr 2 Ru0 4 



The inter-plane coupling in Sr2Ru04 leads to a super- 
conducting ground state driven by the d zx and d yz cor- 
relations. The d zx and d„ * c orrelations are described by 
the model derived in Ref. [32| and applied in Sec. [II. The 
pair order parameter is the expectation value of orbital- 
singlet, spin-triplet Cooper pairs. The gapless fluctu- 
ations of the multiple order parameter components are 
dominantly two dimensional. In this section the physical 
implications of the model are compared with the proper- 
ties of Sr 2 Ru04 determined through experiments. 



and consequently 



F, 



T~T C 



V F Q4l + NV ^P 2 + NbP 4 . (54) 

A* 1 c 



By analogy we adapt this description near the phase-tt 
sition while for T — > the free energy saturates aal; 



F M (T = 0) ~ P 



l-(-Bg)/(8f) 



(55) 



The dimerized spin chain is SU(2) invariant with B 2 ^ 




001 Since Bl 
1631 at T 



2tt is not a singular point in the 
similar behavior for other values of 



can be expected at finite temperature. The com- 



parison of Eqs. (54) and ( p5[ ) yields an effective finite 
temperature interaction parameter of 



A. Order parameter: /iSR and excess tunneling 
current 

Minimizing the Landau expansion of the free energy 
Eq. (pi) via (8F A )/(dP) = gives 



it<i 



9V£s 



cir 



t. 



(59) 



The slope „ r, c | f g « 0.1 has been determined from the 
analysis of the specific heat data in Sec. IV B. For T — > 



or t — > 1 Nernst's theorem requires the order parameter 
to saturate. A corresponding fit function to the experi- 
mental data as shown in Fig. |l| yields using Eq. (|sT 



^(T~T c )«47r, 



(56) 



A(T) = 5 - v/l + llT 2 /K 2 K 



(60) 
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FIG. 1. Order parameter as a function of temperature. 
Circles are the excess current form Ref. (axes scaled), 
squares are fiSR relaxation rates from Ref. VLawith subtracted 
normal phase relaxation rate. Full squares are for Muon spins 
polarized in the x-y plane, open squares are scaled with 2/3 
for Muon spins polarized perpendicular to the x-y plane. The 
phenomenological fit 10~ 2 (5 - \/l + 11T 2 ) [Eq. @] is con- 
sistent with the linear temperature dependence of the order 
parameter predicted in Eq. 



and is consistent with Eq. (|g) for T c = 1.48 K. Note 
that the determination of T c from experiments for an 
order parameter linear in t is slightly different than in 
the case w ithout gapless fluctuations where A ~ y/t, see 
Sec. EVB . 

The squares in Fig. |l| show the increase of the /zSR re- 
laxation rates AAj_ (open symbols, Muon spin polarized 
perpendicular to the x-y plane) and AAn (full symbols, 
Muon spin polarized parallel to the x-y plane) below the 
superconducting phase transition from Ref. |l^. The in- 
crease of the /zSR relaxation rate is a consequence of the 
gapless magnetic fluctuations induced by the spin-triplet 
Cooper pairs in the superconducting phase and described 
by fi s . The magnetic excitations can scatter dynamically 
off the Muon spins. The resulting coupling to thcJ\luon 



spin must be described in the dynamical limit asE£^ 



A, = Bl t- 1 . 



(61) 



Bq is the field at the Muon site, t^T 1 is the Cooper pair 
spin scattering rate. FipiH-the absence of any change in 
the static susceptibilityti-3'Ej can be concluded that Bq is 
the same in both the normal phase and the supercon- 
ducting phase. The absence of the f2 s fluctuations in the 
normal phase yield A s 3> A„ , where A n is the Muon relax- 
ation rate induced by the "normal" electrons. Since the 
contribution from the Cooper pairs is ~ A S P/P while 
that from the "normal" electrons is ~ A n (l — P/Po) one 
has 



The results from /xSR are consistent with the linear tem- 
perature dependence of the order parameter (Fig. |l|) . 

In the presence of an external magnetic field the mag- 
netic components of the orde r parameter are Zeeman 
split as discussed in Sec. Ill C . The spin scattering rate 
t~ and thus A s are exponentially suppressed. This is 
consistent with the experi ment al observation.113 

As mentioned in Sec. HIE the spin-flavor coupling 



terms in Eq. fl23| ) could in principle lead to a gap in the 
magnetic excitation spectrum. On the other hand such 
a coupling would inevitably add a magnetic component 
to the flavor degrees of frcedpm. Since the SO (2) flavor 
fluctuations remain gaplessEJ the Muon spin would then 
couple dynamically to the magnetic moment of the fla- 
vor mode leaving the physical picture of the origin of the 
increase of the /iSR relaxation rate unaltered. 

From AA|| w 2/3AAj_ it can be concluded that the 
spin-one Cooper pairs are in a slight easy plane configu- 



ration 



B 



0.8S||. Note that AA| 



B\ and vice 



123 From the critical field analysis in Ref. the 



versa 

in-plane moment can be determined more precisely to be 
enhanced by a factor 1.23 with respect to the out-of-plane 
moment. 

The excess current across a normal-superconducting 
point contact is proportional to the superconducting gap 
at le ast in s-wave superconductors.1123 As discussed in 
Sec. II C and Ref. the order parameter in Sr2Ru04 



has extended s-wave character and in Sec. IV E the su 



AA = (A s - An) P/Po 



(62) 



perconducting gap in Sr2Ru04 is determined to be 93% 
isotropic. The full circles in Fig. [l] are excess current 
data from Ref. [l?] and are in striking agreement with the 
linear temperature dependence of the order parameter. 
The data are scaled to match T c and the slope. 

Equation ( |60| ) gives an energy gap of A|t^o ~ 0.33 
meV. The results from the differential resistance versus 
bias voltage measurements across point contact junctions 
depend on the model used for their analysis. A p-wave 
analysis without gapless excitations yields AJw^o ~ 1-1 
meV, a s-wave analysis A|t^o ~ 0.25 meV.EZl Both re- 
sults are close to the value extracted here. 

The differential resistance data have been argued to 
be inconsistent with an isotropic s-waye. superconductor 
in the absence of gap less excitations.OEa The theoretical 
approach in Ref. 104 supports this conclusion but ignores 
contributions from the d zx and d yz orbitals. An analysis 
of the Andreev reflection spectroscopy data in the pres- 
ence of the herein discussed order parameter including 
massless fluctuations is desirable. 

he p-wave approach proposed in the 
3 the order parameter d(k) = z(k x ± ik y ) 
breaks time reversal symmetry and thus-can also explain 
the change in the Muon relaxation rateEjfl at the super- 
conducting phase transition. The resulting temperature 
dependence of the order parameter enters via the local 
field and is consequently y-iVi consistent with the dif- 
ferential resistance analysisO but not with the measured 
linear temperature dependence of the excess current. On 
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FIG. 2. Specific heat over temperature as a function of 
temperature. Symbols are from Ref. ^(j. T cm = 1.48 K is the 
"mid transition" value. The actual transition is at T c = 1.52 
K. In the normal phase the specific heat is linear (All line) 
consistent with Fermi or Luttinger liquid behaviorEn For T 
just below T c a linear increase with t = 1 — T/T c is observed 
consistent with an order parameter linear in t [Eq. (pqjj. 
For T < 0.3 K the specific heat is quadratic in temperature 
[dot-dashed line, Eq. ffity]) consistent with two-dimensional, 
gapless order parameter fluctuations. Dotted line: sketch of 
the BCS contribution. 



is consistent with the measu red superconducting energy 
gap as discussed in Sec. IV A . One finds TJ^/sq « 22 K 2 . 
Note that the "mid transition" temperature T cm — 1.48 
K is slightly lower that the critical temperature T c = 1.52 
K determined here. 

For slightly lower temperatures, i.e., t > 0.05, the neg- 
ative, higher order terms in Eq. ( |63| ) quickly become dom- 
inant. For temperatures 0.7 > t > 0.05 the experimental 
specific heat is given by a superposition of the exponen- 
tial term from the gapped electronic charge, spin, and 
spin flavor channels [Eqs. (4£) and ([50])] and from the 
non-linear sigma model [Eq. (47)]. 

The appropriate sine-Gordon actions for the electronic 
channels are given by Eqs. @, ©, and (|l]). With 
finite Mfj, ^ they yield a BCS-like, exponential specific 
heat contribution as sketched be the dotted line in Fig- 
ure ^| The fluctuation contribution from the non-linear 
sigma model Eq. ( |47j ) to the specific heat in Eq. ( |63| ) is 
negative ~ —t for t — > 0. Consequently the gapless fluc- 
tuations account for the jump in the specific heat at T c 
that is lowec. than anticipated by a gapped mean-field 
calculation.EJ The physical interpretation is that a sig- 
nificant contribution to the entropy is contained in the 
gapless fluctuations of the internal degrees of freedom of 
the order parameter. 

For small temperatures t > 0.8 the contribution from 
the electronic channels is negligiblec£l and the quadratic 
temperature dependence of the non-linear sigma model 
Eq. (M2|) dominates the specific heat. 



the other hand, the excess current might be dominated 
by surface effects not considered here. A more obvious 
possible shortfall jof the p-wave approach are the neutron 
scattering resultstJ that do not show any signature of a 
gap down to energy transfers of u> > 1 meV suggesting 
that |2A < 1 meV which is inconsistent with the gap 
valueu-3 of AIt^o ~ 1.1 meV discussed above. 



B. Specific heat 

The normal phase thermodynamic properties have 
been described satisfactorily in the framework of the 
model outlined in Ref. |32|. The specific heat C^ ot s» 
37.5 ^a^j T is linear in temperature for T c < T < 30 K 
as is shown by the solid line in Fig. |[ 

In the superconducting phase near the phase transi- 
tion, i.e., t = 1 — T/T c -C 1, the free energy is given by 
Fa in Eq. ( |58| ) while the temperature dependence of the 
superconducting pair density P is given by Eq. (^) . The 
resulting specific heat C s = -Td 2 F&/ (dT 2 )\t >0 is 



C. 



s|«l 



C\, 



V tff 



6 Vgsl 



(63) 



For V = 6 K (see Sec. II A) and T c = 1.52 K according to 



the fit in Fig. (dashed line) one has 



3<(2) T 2 



3.4 



(v cS ) 2 



(64) 



Here £(n) is Riemann's zeta function and the prefac- 
tor 3 is the-jnultiplicity of the channels fj, = s, f that 
contribute£23 From the fit in Fig. || (dash-dotted line) 
one extracts v e s m 22 K. The resulting value of so ~ 4.7 
is only an order of magnitude estimate since the accurate 
prefactor in Eq. ( fisj ) is not known. If one assumes the 
magnetic channel to be gapped the numbers are v e g ~ 38 
K and sq w 2.7. 

The value of 2 v e ff w 44 — 76 K must be compared with 
the quasi one-dimensional magnetic excitation velocity 
determined in Ref. |3^ as v c s ~ 10 2 K albeit outside the 
Fermi liquid regime. The discrepancy must be attributed 
to effects of the two-dimensional coupling and the pair 
correlations in the superconducting state. The fact that 
the numbers are of the same order of magnitude is a non- 
trivial consistency check of the theory. 

The specific heat has been reproduced fairly ,. well, 
within the p-wave approach by different groupsEjtHjEJ 
by modeling the gapless excitations via introduced line 
nodes or a two-gap scenario,c3 where one of the gaps is 
very small. Since a number of experimental jyajaes re- 
quire a rather homogeneous in-plane gapE3E3'a'E3 only 



6W^ 



horizontal line nodes are possible. The. latter rrq.uire the 
fine tuning of a number of parameters .l3cZI : E£3l3 The fi- 
0.1 which m te slope near T c has not been addressed yet. It would 
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FIG. 3. Normalized Inverse stiffness p° s / P T(T) (full line) 
from Eq. for pi « v e g. The broken line shows the asymp- 
totic part v cff /pj mp . 



be interesting to repeat those approaches using the pair- 
ing potential in Eq. (^) including the gapless fluctuations 
of the internal degrees of freedom of the order parameter. 



C. Stiffness and effective saddle-point regime 

The stiffness near the phase transition, where the 
amplitude-Goldstone mode coupling is relevant, has been 

determined in Sec. |IIIE2| as pf llp ~ P{X 2 )A V 



T- 



At low temperatures it has be determined within the non- 
linear sigma model Eq. ( fi2j ) as ~ P A„\t^o- Since at 
low temperatures the exact relation between P and the 
gap is not known [c.f. Eq. (ftO| )] a precise numerical de- 
termination of pi is not possible. For sufficiently strong 
coupling the reasonable assumption is that the stiffness 
is of the order of the typical energy scale of the system, 



i.e. 



Ueff (c.f. Ref. 



To obtain a qualitative picture of the temperature de- 
pendence of the stiffness the approximate interpolation 
formula 



P T(T)^j(p^ + ( P r p ) 2 ^v cS \i + 



T*v 2 



off 



(65) 



can be applied. The full line in Fig. || shows the normal- 
ized inverse stiffness from Eq. ([35]) , the broken line is the 
asymptotic part v e s/Ps mp - The temperature dependence 
of A(T) has been modeled via Eq. (|60|). 

At intermediate temperatures 0.8 > t > 0.05 the stiff- 
ness increases as ~ t~ 2 and the integrals in the non- 
linear sigma models need to be cut off as in Eq. ([l7|). For 
t < 0.05 one has v e f[/p^ mp < 10 -3 which is sufficiently 
small to justify the effective saddle point approach to 
the fluctuation action Eq. (48) and the resulting Landau 



expansion of the free energy Eq. ( p8[ ) . 

For t > 0.05 the effective saddle point approximation 
Eq. ( |4"g| ) is not applicable any more. The third order 
term in Eq. (|58| ) then is not present in that form any- 
more which enhances the importance of the higher order 



terms in Eq. (63) and contributes to the rather sharp 
maximum of the specific heat. Note that even in BCS 
theory the specific heat is much more sensitive to the 
higher order terms that stem from the opening of the 
single particle gap than the temperature dependence of 
the order parameter itself. 



D. Knight shift, static susceptibility, NQR, and 
universal thermal transport 

The explanation-©f the observed absence of a change 
in the Knight shifttJ and the magnetic susceptibility^ in 
the superconducting phase is trivial within the present 
model. All order parameter components carry spin one 
and thus no electronic magnetic moments are lost in the 
superconducting phase. 

The relaxation times determined by NQR are consis- 
tent with gapless, two-dimensional fluctuations in the su- 
perconducting stated The qualitative agreement with the 
model presented here is obvious. 

The thermal transport also is consistent with the pres- 
ence of gapless fluctuations below T C M The observed 
universal linear temperature dependence for T — > is 
determined by the long wavelength excitations of the in- 
ternal degees of freedom of the order parameter and is 
consequently determined by the properties of collective 
state much rather then by the number of impurities in 
the system. 



E. Order parameter symmetry 

Equation ( |3l"| ) reveals the spatial order parameter 
anisotropy 



(P) « (Sl t ) 



(sin( 1 /7r<fr 
(sin(y / 7T^f - 



V2x)) 



(66) 



As discussed in Sec. [II C the pair correlations for (Of x ) 
are larger along x or [110] and for (Of,y) are larger along 
y or [TlO]. The simplest parameterization is 



(o f ) = 



(o f 



V2 



sin 2rj 



(O 



r.i 



y/2 \-l 



1 



(67) 



where the angle r\ is measured with respect to the [100] 
axis. Each component has a two-fold symmetry while the 
resulting spatial symmetry of the total order parameter 
reflects the four-fold axis of the underlying lattice. 



O f )| = VRo) 2 + (Of,i) 2 sin^tj 



(68) 



Oriented magnetic fields and tempe ratur e gradients 
break the SO(2) symmetry of Of (Sec. IIIC) and reveal 
the two-fold symmetry oi-the components. Thermal con- 
ductivity measurements 2 .!! suggest that (Of i i)/(Of i o) ~ 
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0.02. A more careful analysis of the upper critical fieldsE^I 
in Ref. H yields (n f) i)/(O f)0 ) = 0.034 which leads to a 
spatial anisotropy of the order parameter components of 
7%. The the four- fold symmetry of the total order pa- 
rameter is then 



Kn f >|«<fif, > i 



(«f,0> 



sin 2r/ 



(69) 



which is also consistent with the experimental 
observations M 

The experimentally implied existence of two order pa« 
rameter components with a slight spatial anisotropyO 
thus follows quite naturally out of the degenerate saddle 
points in the flavor channel discussed in Sec. IIIB. 



V. RELATION TO THE p-WAVE APPROACH 

A large number of experiments in superconducting 
samples of Sr2Ru04 have been interpreted in terms if the 
time reversal symmetry breaking p-wave state with_S¥m- 
metry E u and order parameter d(k) = z(k x ±ik y )E3'U As 
has been discussed throughout the manuscript most ex- 
perimental results find alternative interpretations within 
the present approach. Examples are 

(i) the obseriied— density of states below the super- 
conducting gapEaLilB that has persistently been inter- 
preted as an indication for line-nodes in the gap func- 
tion while it may well originate from any kind of gapless, 
two-dimens ional excit ations with linear dispersion (Sees. 
EITe| , [IV B| , and [TV Eft ; 

(ii) the c-axis pressure dependence of T c that was in- 
terpreted in Ref. ^ as a consequence of a shift of the 
7 band closer to the van Hove singularities near the M 
point while in Sec. || it is viewed as the simple enhance- 
ment of the inter-plane coupling; 

(iii) the strong in-plane coupling of strain to the super- 
conductivity which was interpreted as a possible indica- 
tion of an oder parameter of the form d(k) = z(k x ±ik y ) 
while the anisotropy of the components of J~2f can account 
for the same effect for symmetry reasons (Sec. [II E[ ) ; 

(iv) the change in the /xSR relaxation rate across the 
superconducting phase transition that was interpreted as 
an indication, for a time-reversal symmetry breaking pair- 
ing stateE3S while it may as well stem from the scatter- 
ing of gapless magnetic excitations off the Muon spin 
(Sec. [EVA| ); 

(v) the sensitivity of the superconductivity to non- 
magnetic impurities which has been interpreted as 
an iridic atiipa, for higher angular momentum pairing 



pairing state similar to the p-wave E u state. The fla- 
vor expectation values (Of) in Eq. ( |3l| ) possibly con- 
tain contributions that are odd in real space, i.e., 
(sin y/2y cos y/nfa) ^ and (sin y/2x cos y/n<j)f) ^ 0, and 
consequently a representation 



{ftp} OC 



(sin y/2y cos i/tt(/){) + i(sin y/2x cos y/nfa} 
(sin y/2y cos y/nfa) — i(sin y/2x cos y/nfa) 



(70) 



ID 

states 

metry of the orbital-singlet pairing state as discussed in 
Ref. |J. 

These similarities between the approaches appear not 
very surprising when considering that the spin-flavor 
phase space discussed herein includes a time-reversal 



can be constructed with fi p G {f2f <£> Jl s }. The magnetic 
configuration of tt p depends on the details of the spin- 
flavor coupling te rms a nd an easy-plane configuration as 
discussed in Sec. IV A and in Ref. |4l] is quite conceiv- 
able. The present approach with mixed-orbital, singlets 
can thus be generalized to include the p-waveBO E u 
picture. 

The conceptual advantage of the present approach is 
its derivation from a microscopic model. It allows to 
study interaction effects like the presence of two degen- 
erate saddle points in detail — at least on a qualitative 
basis. Phenomenologically the important difference is 
that here Hund's rule coupling is emphasized over effects 
from spin-orbit coupling as discussed in Sec. II A which 
allows for t he de generacy in the internal degrees of free- 
dom (Sec. HIE). This leads to a coherent description 
of a multitude of observations in Sr2RuC>4 including the 
experiments discussed herein as well as the normal phase 
magnetism (Ref. [52]) , the interplay of the lattice symme- 
try, impurities and interaction effects (Ref. 34), and the 
magnetic field dependence (Ref. |4l]) . 

Resulting specific advantages over the p-wave E u de- 
scription are 

(i) no need to construct horizontal line nodes or two- 
gap states which require parameter fine tuning and which 
have not been confirmed experimentally; 

(ii) no inconsistency between the gap value extracted 
from Andreev reflection spectroscopy with respect to 
neutron scattering results (Sec. IV A); 

(iii) the correct two-component order parameter sym- 
metry as determined via therm al tra nsport and upper 
critical field measurements (Sec. IV E); 

(iv) the description of the t emperature slope of the 
specific heat near T c (Sec. 1VB) as well as the magnetic 
field dependence of the specific heat near H c (Ref. f4l| ). 

To obtain a more conclusive discrimination of the state 
promoted in this series of papers and the p-wave E u 
state we pxqpose to study a detailed Landau-Ginzburg 
functiona l 10 ^ that contains in its parameter space the 
degenerate state described herein as well as the special 
caseL30 given in Eq. ( |70|) . Experimentally, the presence 



while it follows directly out of the sym- f the fl { mode might be detectable via microwaves. 



VI. CONCLUSIONS 

From the approach discussed in the present paper the 
essential physical properties of the superconducting state 
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in Sr2Ru04 can be summarized as follows. 

The body centered tetragonal structure gives rise to 
inter-plane pair correlations in the d yz and d zx orbitals 
enhanced by umklapp scattering processes. The super- 
conducting order parameter is found to be a mixed- 
orbital-singlet spin-one-triplet. The bosonized descrip- 
tion of the in-plane electron correlations is consistent 
with the order parameter and the slight easy plane con- 
figuration of the magnetic moments. The model reveals 
that the oder parameter has two slightly anisotropic spa- 
tial components. 

The absence of a change in the Knight shift and the 
magnetic susceptibility in the superconducting phase fol- 
lows trivially from the spin-one Cooper pairs. The ex- 
perimental thermal conductivity and upper critical fields 
are consistent with a spatial anisotropy of the order pa- 
rameter components of 7%. 

The different components of the order parameter in 
the spin and flavor sector give rise to two-dimensional 
gapless fluctuations. They are modeled by a 2+1 dimen- 
sional non-linear sigma model. Near the phase transition 
the two-dimensional gapless fluctuations account for the 
linear dependence of the pair density and the specific 
heat on the reduced temperature. At low temperatures 
they yield the quadratic temperature dependence of the 
specific heat, the cubed temperature dependence of the 
NQR relaxation time, and the universal linear tempera- 
ture dependence of the thermal transport. 

The finite pair density of the electrons in the d xy 
orbitals is induced by the inter-band proximity effect. 
There is only one phase transition in the absence of mag- 
netic fields and one single particle gap. 



relevant orbital-singlet spin-triplet channel are given as 



P s t = — sin y/n^st flp <B> O s <8> Of 



(71) 



Recall that the flavor degrees Of of freedom stem from 
the symmetry of the bosonized action which suggest a 
two- fold degenerate saddle point. 

The superconducting state that is discussed in this pa- 
per involves the breaking of various continuous symme- 
tries. The broken gauge symmetry is parameterized via 
the charge phase (O p ), the SO (3) symmetry of the mag- 
netic moment of the spin-triplet Cooper pairs enters via 
(fi s ), and the broken SO (2) symmetry of the flavor chan- 
nel is given by (Of). The expectation value of the order 
parameter is given in vectorial representation via 



(p s t ) = p o p <g> o s ® o f , 



(72) 



can be absorbed in O p via a 



where P is the order parameter amplitude as defined in 
Eq. fl34| ) and O p (r) are normalized according to Eq. (|3~" 
The phase of (s'm^/TT(j) s { 
gauge transformation. 

The ground state is degenerate with respect to ro- 
tations of n M giving rise to long wavelength gapless 
excitations — the Goldstone modes. In order to describe 
the dynamics of the system in the ordered phase in proves 
useful to introduce pair operators that are related to 
those in Eq. O by a shift of the pair expectation value, 
i.e., P = P s t ^(Pf). The mean-field decoupling Eq. @ 
then can be rewritten as the exact relation 



p V ps = p V{ ps 



P s t (Pt )*-\(pi )l 2 + p'p 



(73) 
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APPENDIX: DERIVATIONS FOR THE ORDER 
PARAMETER FLUCTUATIONS 

The expressions that are used in the manuscript to de- 
scribe the fluctuations of the internal degrees of the order 
parameter can bp dsijiijed within standard perturbative 
approaches. oLjOEj'Ej Withi n the vectorial representa- 
tions introduced in Sec. Ill B the pair operators in the 



Using Eq. ( |73[ ) the orbital-singlet spin-triplet part of the 
pair Hamiltonian Eq. (^) decouples into two parts H v = 
H m f + Ha, where the mean-field part is given by Eq. ( |l6| ) 
or, in its bosonized version, by Eq. (|35|). The fluctuation 
part is 



N 



p (g). 



(74) 



The partition function of the whole system can now 
be expxessecL ._. ida a standard perturbative 

a pprna cW 10 tNj 8 4Mf 5 l as 



z 



" J V[X,X*} exp{-S R [X,X*}}. 



(75) 



The mean-field action is composed of Eqs. (|39|), (p8|), and 
([ll]) as S m { = S° s + S°f + S c and accounts for the physics 
contained in the saddle point including the mean-field 
temperature dependence of the gap and the Goldstone 
modes as discussed in Appendix A below. The ther- 
modynamic averaging in the amplitude fluctuation part 
~ S$[X, X*] is performed with respect to the mean-field 
system Smf- The vectorial character of th*j-»air opera- 
tors leads to Hubbard-Stratonovich ficldsBEJ'c2l X and 



1G 



X* which also have vectorial character and can be pa- 
rameterized as 



X = X fi„ ® fL ® f2 



r ■ 



(76) 



The integral measure contains the averaging over 
the aHji^liaij components of the Hubbard-Stratonovich 
field/" 



V[X,X*} = Y[ Jp\v q 



' dX q n dfl p dfl s dflf 



(77) 



2n An 2n 



The to n are bosonic Matsubara frequencies. The av- 
eraging over the angular components corresponds to 
an aver.aging over the degenerate ground states of the 
system.EJ'E-3 



A. Goldstone modes 

In the present approach the Goldstone modes that cor- 
respond to the broken local spin rotational, flavor, and 
gauge invariance are implicitly contained in the mean- 



field system defined by the actions S m f = S[ 



off 



Sf + S c 



This becomes obvious when considering that the mass 
terms in Eqs. (|39|), Q3q), and pl| ) are decreased when 
the operator fl^ is not parallel to the mean direction 
ftp. A direct derivation of the non- linear sigma model 
in Eq. ( pl2| ) that describes the dynamics of the Goldstone 
modes within the present approach would require non- 
Abelian bosonization in order to maintain the symmetry 
of the Hamiltonian explicitly.OE3 In the light of the well 
established results in the literature the involved explicit 
derivation herein appears obsolete. 

A simple qualitative argument that makes the un- 
derlying physics transparent to the reader shall be 
discussed instead. The finite local angle ?7 M (r,r) = 
Z[fi M (r, t), fi M (r, t)] increases the energyEj of the system 
since the mass term determining the gap is lowered. Con- 
sequently there is a restoring force which can be assumed 
linear in 77^(7-, r) if the angle is sufficiently small. Small 
fluctuations of ^(t^t) are consequently expected to be 
harmonic. Their spectrum is gapless since the restoring 
force vanishes for r? M (r,r) ^£L.-T-he excitation velocity 
should be of the same orderEJ'ESIiiil as the velocity of the 
gapped amplitude mode (Appendix B). Finally, the local 
constraint of I^^t)! 2 — 1 distinguishes the angular 
fluctuations from other harmonic systems and leads to 
the effective description via th e non-linear sigma models 
given by Eq. @ in Sec. [HIE. 



At sufficiently small temperatures and energy scales 
T, u> -C A the stiffness of the non-linear sigma 
model describing the angular fluctuations has been 
determinedOO to scale with the square of the order pa- 
rameter amplitude ~ P 2 , albeit with correctionsBIl! that 
depend on the specific correlations of the system. As 
will be discussed in Appendix B2-the gap in the ordered 
phase plays the role of a cutoff.O This cutoff becomes 



small as the gap becomes small and has important phys- 
ical implications.!!^ 



B. Gaussian amplitude fluctuations 



The integral term in Eq. (|75| ) contains the fluctuations 
P of the system involving changes in amplitude of the 
order parameter. The amplitude fluctuation part is ex- 
panded up to second orderHBa in the complex fields X 
and X*. Considering that per definition (P) = one 
has 



S fl [X,JT] - 



V q \x { 2 P i(q,iu n ) \X(q,iuj n 



(78) 



For the relevant umklapp scattering near q ~ the po- 
tential has been set V q w — \V q \. The pair correlation 

(p) 

function \2 ± needs to be calculated with respect to S m {. 

Since S m { describes a gapped system in the ordered phase 

(p) i (P) i 

X2.±(°,0)| T<T < X2,±(°)°)|r=T anc * tne action in Ec l- 
@ is stable. ° 

Limiting the description to the low energy long wave- 
length properties one may expand 



Wx<f2(g,w„)«|Vo| 2 x£ P M°) 



2,± v 



A x ql + A v q, 



A x 
(av x ) 



A*. 



{avy) 2 



(^„) 2 , (79) 



where A u and v v arc defined after analytical continua- 
tion iui n «-> \\t[i c ^q{lu + ie) in Eqs. (43| ) and (fl4|), respec- 
tively. As discussed in Sec. HIE 1 one can assume for 



the low temperature and low energy regime the system 
to be isotropic in the x-y plane with A v = A x = A y and 

v u =v x = Vy. 

Using Eq. ((79|) and transforming Eq. (|7S| ) to Euclidean 
space representation the action becomes 





d'r I dr \Vo\ |1 



b|[l 



1^01X^2(0,0) 



d 2 r 



dr 



E 

v=x,y 



A„ 



\3 T X\ 



\d v X\ 



(80) 



The q dependence of the term linear in \V q \ has been 
neglected. Note that in Euclidean space \X\ 2 = X 2 
because of the local constraints |i~2 M (r,T)| 2 = 1. An- 
other important consequence of the local constraint is 
that J7 M i9 T r2 M = Clpd^flp = 0. Consequently angular 
and amplitude contributions to the derivatives in Eq. ( |80| ) 
decouple for v = x, y, r as 



\d u x\ 2 = (d u xy 



x 2 J2 ( d - n 

fl—p,S,f 



%) 2 



(81) 
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Now the different modes around the saddle point can be 
discussed in detail. 



1. Amplitude mode 

The part of Eq. ( |80| ) that describes the amplitude ex- 
citations is in Fourier space 



Samp[^, X*\ — (3 



M(l-|V |4f2(0,0) 



\X(q,iu n )\ 2 
(82) 



and describes Gaussian fluctuations around the static 
saddle point. 

In order for the action Eq. ([82]) to be consistent with 
the mean-field pair correlation function estimated by Eq. 
(fl5|) one must require that 



P{X(q,iw n )X(-q,-iu) n )) = XBCs(l> iu; ri) 



(83) 



After analytical continuation iuj n <-> lim e _>o( w + * e ) an d 
comparing Eqs. (|82]), (|S3"1), and Xbcs m one can 

identify 



-4, 



A/ 



as well as 



\V \(l-\V \ X ({l(0,0))*j o 



(84) 
(85) 

(86) 



The q dependence of /(g) s» /(0) = /o has been ne- 
glected. The dispersion of the amplitude mode then can 
readily be read off from Eq. (|82j) as 



^amp 



v eS Ja 2 q 2 - A 2 /v 2 cB 



(87) 



and has the desired shape with an excitation gap A un- 
derlining the self-consistency of the approach. 

The identical self-consistent results are obtained by di- 
rectly identifying 



Wo 



l-\V q \xZ2(q,iu n ) 



Xbcs(<?>* w «) ( 88 ) 



in the action Eq. ( |T8|) which again shows the stability of 
the latter since Xbcs(0' 0) > 0. 



2. Amplitude- Goldstone mode coupling 

Starting from the fluctuation action Eq. ( |80| ) , including 
the identity Eq. J8l|), using the identifications Eqs. (H) 
and fl85|), and transforming to Fourier space the action 

s^[n„n;]=pJ2 E A » E * 2 (^«) 

q,u n u=x,y (j.=p,s,f 



x (aq„) 2 + 



\n^(q,iu) n )\ 2 (89) 



is obtained for the angular modes. Except for the pref- 
actors it is equivalent to the non-linear sigma model Eq. 
( pl2| ) used to describe the low temperature, low energy 
excitati ons o f the angular Goldstone modes [Appendix 
A, Sec. |IIIE| . 

The action Eq. (|8^) is the lowest order amplitude- 
Goldstone mode coupling term. This becomes evident 
when replacing the amplitude of the pair fluctuation field 
X by its expectation value, i.e., X — > (X) — 0: in the 
absence of amplitude fluctuations Eq. (|89| ) does not con- 
tribute. This is reminiscent of the decoupling of the pair 
operators into angular and amplitude parts through Eq. 
( [73| ) where the angular fluctuations that give rise to the 
Goldstone modes have been included in the mean-field 
action S m { [Appendix A]. 

In a rigorous approach the amplitude fluctuations have 
to be integrated outEj in order to determine the contri- 
bution of the action Eq. ( |89| ) to the free energy of the sys- 
tem. In the light of the local constraint |r2 M (r,r)| 2 = 1 
the treatment of the resulting terms is not evident. In- 
stead, on the mean-field level one can estimate the im- 
pact of the amplitude- Goldstone mode coupling by re- 
placing the square of the amplitude fields by their ex- 
pectation value, i.e., (3X 2 (q,ioj n ) — ► (3(X 2 (q,iuj n )) = 

XBGs(9,^n). in E q- ©■ 

Xbcs i s gi ven in Eq. ( f45| ) and diverges for uj 2 — v 2 q 2 — » 
A 2 . Consequently the excitation gap A indeed-Miarks the 
ultraviolet cutoff discussed in the literature. EiiE2l Given 
the mean-field relation between the excitation gap and 
the coherence length A/v e s ps £bcs the cutoff translates 
into a short range or large wavevector cutoff a|oL< A/v v 
in strict analogy to the discussion by SchulztS for the 
magnetic order in the repulsive Hubbard model. 

For sufficiently low temperatures T <C A the thermal 
occupation of states near the cutoff energy can be ne- 
glected and the gapless excitations are well described via 
the n on-linear sigma models given by Eq. (^]) in Sec. 
IIIE . Since in this case the amplitude mode has a large 
gap compared with temperature, contributions from the 
amplitude mode and the amplitude-Goldstone mode cou- 
pling term are negligible. 

Closer to the phase transition corrections from Gaus- 
sian fluctuations become relevantES and the energy and 
wavevector cutoff will play a crucial role and needs to be 
included in the action Eq. (|8^). In continuum representa- 
tion for the momenta this leads to Eq. (^7j) in Sec. Ill E 2 . 
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Note that for small energies and wavevectors the us and 

q dependence of XBCsil^^n) ~ XbcU * ) = /o/ A and 
\fl^(q,iuj n )\ 2 ~ |^| 2 is not expected to alter the phys- 
ical result and consequently have been neglected in Eq. 
©■ 



3. Limits of the applicability 

It is well known@ that in the critical region around the 
phase transition the Gaussian approximation for fluctu- 
ation corrections as derived in Appendix B breaks down. 
The contributions are divergent and higher oder terms 
have to be included in a non-perturbative fashion. The 
critical region can be estimated by the rounding of the 
measu red sp ecific heat near the phase transition (Fig. ^ 
in Sec. IV B ) and yields a temperature interval of T c ±0.02 
K. The results herein including Eq. (|5^) are only appli- 
cable for T < T c — 0.02 K or for reduced temperatures 
\t\ > 0.01. Note that within the framwork of the third 



order phase transition considered in Sec. [VB the deter- 
mination of T c and thus of the critical region differs from 
the usual "mid transition .point" estimate. 

The possible impactE3'E3H_^of additional topological 
terms is addressed in Sec. Ill Ej . 
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